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Comment on "Two-spinon and four- 
spinon continuum in a frustrated ferro- 
magnetic spin-1/2 chain" 

Recently Enderle et al. reported an inelastic neutron 
scattering (INS) study of the dynamic spin susceptibil- 
ity Imx(w,fc) for LiCuV0 4 Q. Therein they claim that 
(i) LiCuVC>4 is well described by two interpenetrating, 
weakly ferromagnetically (FM) coupled Heisenberg an- 
tiferromagnetic spin-1/2 chains (HAF), (ii) the obtained 
exchange integrals J; (NN and NNN inchain couplings 
Ji=— 1.6 mcV, J2=3.56 mcV, diagonal interchain cou- 
pling J5=— 0.4 meV in the (ab)-planc, a=— J2/J1 ~ 2.2) 
agree with those from an analysis based on spin-wave 
theory (SWT) 0], ( iii) the observed INS intensity above 
10 meV belongs to a 4-spinon continuum (4SC). Apply- 
ing exact diagonalization (ED) and DMRG methods to 
fit INS and magnetization M(H) data, supported by in- 
dependent microscopic methods Q, we will show that 
the claims of Ref. [l| are not justified and that LiCuVC>4 
exhibits a <1, i.e. strong coupling of the HAF, at odds 
with (i). For possible spin nematics and Bose condensa- 
tion of 2-magnon bound states in LiCuV04 0, [B| precise 
knowledge of the coupling regime is of key importance. 

Starting from a 2D model, the authors suggest an ef- 
fective ID-model with J eff ,i=Ji + 2Js =— 2.4 mcV and 
a renormalized J e£ f > 2=2 J2/71", i.e. a~1.4. But the ap- 
plied perturbative method is designed for a 3>1. Fur- 
ther problems occur for the dispersion of spin excitations 
ui(k) (dark red curve in Fig. 2 of Ref. 1). Near fc=7r/4 
the local maximum ft gives 4.84 meV. For a=1.42 and 
J ei s,i=— 2.4 meV one has 51=4.36 (4.79) mcV, only, ac- 
cording to our dynamical DMRG (ED) calculation with 
L=96 (28) sites of hnx(u,q). Fitting our DMRG results 
for 0.5 < a < 2, the general constraint for fl(a) reads: 

4.84meV = = J cfl .i (0.573 - 1.702a + 0.0109a 2 ) . (1) 

More strikingly, fitting the M(£f)-data at T=1.6 K [10 
by our DMRG (T=0, L=512 sites), we found a =0.75 
(see Fig. la). We obtain ff„=41.6 T, ,g=2.27 and in ac- 
cord with Ref. 5 d^/d-g^O.39 at low H, where H S (M B ) 
is the saturation field (magnetization). Then, Eq. (1) 
gives J off ,i=— 6.95 meV and J effj 2=5.2 meV. The ID-set 
of Ref. 1 yields a too small i? B =37.3 T at g=2.27 or a too 
small 5=2.03 for H S =A1.6 T, only, and clearly too high 
M above 0.6H B (see Fig. la). If the SWT-fit is meaning- 
ful, Ji is strongly renormalized but J e g,2 is close to its 
bare value of 5.6 meV both at odds with Refs. 1,2. 

With our fitted J G «2-value, almost the whole region 
ascribed to the 4SC [l[ is covered now by the extended 
2-spinon continuum (2SC). The 4SC should be looked for 
at k = 0.5 above 16.3 meV (see Fig. lc) i.e. in a region 
which has not been measured yet. Since most of the INS 
intensity below 16 meV belongs to the 2SC, the size of the 
4SC enhancement compared to that of a HAF is not yet 
settled and the 4SC in Fig. 4 of Ref. 1 is overestimated. 
Figs, lb, lc show that the dispersion of the INS peaks is 



insufficient to find a unique J1-J2 set. INS intensities or 
other quantities like M(H) must be analyzed, too. Our 
set explains the larger INS intensity above 9.5 meV as 
compared with that in Ref. 1 (see the boxes in Fig. 1). 

To conclude, a weakly coupled HAF model with small 
Ji| < 5 mcV as in Ref. 1 is not justified for LiCuV04 
whereas a strong coupling regime with |Ji|></2 where 
both J's exceed significantly those of Ref. 1 is consis- 
tent with the INS, the magnetization, and realistic micro- 
scopic models of edge-shared cuprates For strongly 
coupled HAF's, the 2SC is extended to higher u with less 
spectral weight left for the 4SC. Future studies are highly 
desired to refine the main J's and the 4SC as well. 
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FIG. 1. Magnetization of LiCuV0 4 for the ID J-set of Ref. 1 
and our fit compared with experiment (H || c ) (a). Imx(w, k) 
from exact diagonalizations for both sets for a chain with 
L—28 sites broadened with 0.05 meV. Red curves: dispersion 
of spin excitations from main INS peaks 0(b,c). Notice a 
slight downshift near fc=7r/4 of 0.426 meV and 0.676 meV for 
Fig. lb(lc), respectively, within the DMRG (see Eq. (1)). 
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